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1.  INTRODUCTION 


This  manual  is  divided  into  two  parts.  Section  2  constitutes  the  user 
section,  describing  the  overall  flow  of  the  codes,  input  data  requirements, 
and  parameter  relations.  At  the  end  of  the  section,  a  series  of  sample 
problems  is  given  to  illustrate  the  use  of  the  codes.  Sample  problem  1 
exercises  all  parts  of  the  computer  algorithm  and  serves  as  a  check  case  for 
the  installation  of  the  codes  on  a  given  machine.  Section  3  provides  a 
detailed  description  of  all  the  codes  and  subroutines.  The  user  need  not 
make  reference  to  Section  3  or  to  the  analysis  in  Volume  I  of  this  report  to 
successfully  install  and  execute  the  algorithm  on  a  particular  computer 
system. 

A  brief  synopsis  of  the  main  program  features  and  options  are  listed 
below: 

Flexible  body  geometry  -  The  algorithm  can  treat  a  finite-length  BOT  of  any 
asymmetric  cross  section.  The  analysis  assumes  that  the  body  is  open  on  the 
ends.  Certain  degenerate  forms  of  the  BOT  configuration  have  also  been 
analyzed  with  this  algorithm  such  as  flat  plates,  parabolic  and  square 
cylinders . 

Arbitrary  antenna  placement  and  excitation  -  The  algorithm  treats  single  or 
multiple  rectangular  aperture  (slot)  antennas  embedded  in  the  BOT  surface. 
Location  of  the  antennas  can  be  anywhere  on  the  BOT.  Spacing  between  adja¬ 
cent  apertures  can  be  electrically  small.  All  apertures  can  be  asymmetric. 

Surface  currents  -  The  algorithm  outputs  the  surface  currents  on  the  BOT  sur¬ 
face  in  spatial  and  modal  form.  Both  the  magnitude  and  the  phase  of  the  t- 
and  z-directed  components  are  given. 

Choice  of  polarization  and  radiation  planes  -  The  radiated  fields  and  the 
power  gain,  normalized  to  an  isotropic  radiator,  are  given  in  the  user- 
specified  planes  for  6  and  <j>  polarization. 


Arbitrary  choice  of  sampling  points  for  near  fields  -  The  electric  and 
magnetic  fields  (six  components  in  all)  resulting  from  currents  Induced  by 
antennas  on  the  BOT  or  by  incident  fields  are  computed  at  user-specified 
points  in  the  vicinity  of  the  BOT  surface. 

Calculation  of  aperture-coupled  fields  -  All  the  field  components  for  the 
electromagnetic  waves  penetrating  a  passive  rectangular  aperture  on  the  BOT 
are  computed  at  user-specified  sampling  points.  The  size  and  location  of 
the  aperture  are  also  user  specified. 

The  foregoing  features  of  the  algorithm  are  demonstrated  in  the  sample 
problems  of  Section  2.6.  These  features  have  been  tested,  and  the  corre¬ 
sponding  results  are  in  good  agreement  with  available  data.  As  in  all 
modeling,  caution  must  be  exercised  in  applying  this  analysis  to  certain 
problems  such  as  the  computation  of  the  edge  diffraction  from  the  BOT  ends 
and  the  fields  coupled  through  small  apertures.  To  achieve  sufficient 
accuracy,  a  large  number  of  modes  may  be  necessary  in  these  problems. 


2.  USER  SECTION 


The  MM/BOT  computer  codes  are  written  in  FORTRAN  IV,  consisting  of  five 
main-line  programs  (i.e.,  BOTZSS,  BOTINV,  BOTRA,  BOTSCM,  and  BOTSCB)  and  one 
small  utility  program  (BOTSEG)  as  shown  in  Figure  1.  The  five  main-line 
programs  use  the  same  coordinate  data  file  with  minor  deletions,  depending 
on  which  of  the  above  programs  are  run.  In  addition,  two  binary  disk  files, 
ZBOT  and  YBOT,  are  generated  by  the  programs  BOTZSS  and  BOTINV.  All  the 
programs  have  user-oriented  inputs,  which  can  be  easily  generated  for  a  given 
problem.  Examples  demonstrating  the  input  data  requirements  are  discussed  in 
Section  2.6. 

2.1  Implementation  of  Computer  Codes 

All  of  the  BOT  computer  codes  are  written  in  USA  Standard  FORTRAN  IV, 
with  the  exception  of  end-of-file  (EOF)  checks  in  programs  BOTINV,  BOTRA, 
BOTSCM,  and  BOTSCB.  The  EOF  checks  given  in  the  listings  in  Appendix  C  are 
specific  to  the  compiler  used  in  the  program  development  (i.e.,  CDC  CYBER  175 
system) .  The  functioning  of  the  EOF  checks  in  the  present  listings  is  as 
follows: 

IF  (EOF(u) )  a,b 

u  -  unit  number  to  check. 

a  -  statement  label  to  branch  to  if  an  EOF  is  encountered, 

b  -  statement  label  to  branch  to  if  an  EOF  is  not  encountered. 

This  EOF  check  occurs  in  the  listings  (Appendix  C)  at  the  following 

places: 

BOTINV  line  numbers  480,  1310,  and  1870 

BOTRA  line  number  1500 

BOTSCM  line  number  1300 

BOTSCB  line  numbers  1330,  2125,  8310,  and  8970 

The  device  unit  numbering  convention  is  as  follows: 
unit  number  5  -  card  reader 
unit  number  6  -  line  printer 

unit  number  1  -  binary  input  disk  file  containing  alternating 

2 

records  of  lengths  2  and  2(NP-3)  words,  respec¬ 
tively.  NP  is  the  number  of  data  points  describing 
the  BOT  generating  curve  discussed  in  Section  2.2. 


Program  block 


Program  function 


Segments  BOT  generating  curve  (optional) 


BOT  input  data  file  for  the  remaining  programs 


Computes  Zm  n  submatrices  for 
0  <  m  <  NMODE  -  1  and  -  m  <  n  <  m 
(Section  4,  Volume  I) 


Partial  Zf)OT  matrix,  stored  by  submatrices 


Fills  the  ZBOT  matrix  using  mode-dependent 
symmetries  and  inverts  ZboT 
(Section  4.5,  Volume  I) 


Inverted  ZbOT  matrix,  stored  by  submatrices 


Computes  far-  and  near-field  radiation 
(Sections  5  and  6,  Volume  I) 


Computes  monostatic  RCS 
(Section  S,  Volume  I) 


Computes  bistatic  RCS,  aperture  penetration, 
and  near-field  analysis. 

(Sections  5, 6,  and  7,  Volume  1) 

aere-ewi-as 


Figure  1.  Overall  flow  diagram  for  MM/BOT  algorithm. 


unit  number  2  -  binary  output  disk  file  containing  alternating 

2 

records  of  lengths  2  and  2(NP-3)  words,  respec¬ 
tively. 


2.2  BOTSEG  Description 

BOTSEG  is  a  utility  program,  which  can  be  used  to  segment  the  BOT  gene¬ 
rating  curve  using  a  limited  set  of  data  points  (see  coordinate  geometry  in 
Figure  2).  The  required  input  data  and  formats  are  described  below. 


formed  by  a  wedge  and  arc  segment. 

READ(5,1)NPTS,NP 

1  FORMAT (213) 

NPTS  -  Number  of  input  data  points  used  to  describe  the  BOT  gene¬ 
rating  curve. 

NP  -  Number  of  equally  spaced  BOT  generating  curve  data  points  to 
be  calculated.  (The  calculated  data  are  used  as  input  to  the 
remaining  BOT  programs.) 

DO  100  I  -  1,  NPTS 

100  READ(5,2)XTAB(I) ,YTAB(I) ,XC(I) ,YC(I) 

2  FORMAT (4E 10. 4) 

XTAB(I)  -  x  coordinate  of  the  I-th  point  on  the  input  curve 
(meters) . 

YTAB(I)  -  y  coordinate  of  the  I-th  point  on  the  input  curve 
(meters) , 

XC(I) ,YC(I)  -  (meters)  Indicate  whether  the  points  [XTAB(I) , 

YTAB(I) ]  and  [XTAB(I+1) ,YTAB(I+1) ]  are  connected  by  a 
straight-line  segment  or  an  arc  with  changing  radius. 

If  XC(I)  ■  YC(I)  ■  0,  the  segment  is  straight.  Other¬ 
wise,  [XC(I) ,YC(I) ]  is  assumed  to  be  the  center  of  an  arc 
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with  the  above  end-points ,  where  the  radius  changes 
linearly  with  angle  from  [XTAB(I) ,YTAB(I) ]  to  [XTAB(I+1), 
YTAB(I+1) ] ,  subtending  the  angle  of  the  triangle  formed 
by  the  three  points.  Mote,  if  the  three  points  are 
colinear,  the  direction  of  the  arc  will  be  clockwise. 

For  example,  consider  a  cylinder  with  a  cross-section  given  in  Figure  2. 
The  generating  curve  for  this  body  can  be  represented  with  three  points 
(NPTS  ■  3)  as  follows: 

0  0  0  0 

3  13  0 

4  0  0  0 

BOTSEG  is  currently  dimensioned  to  handle  100  points  on  the  input  curve 
and  83  points  on  the  output  curve.  The  NP  data  points  used  to  represent  the 
BOT  generating  curve  determines  the  location  and  number  of  triangle  functions 
and  its  derivatives  (denoted  by  arrays  T  and  TP  in  the  codes) ,  used  to  discre¬ 
tize  the  unknown  currents  on  the  BOT  surface.  As  an  example.  Figure  2b  shows 
a  portion  of  a  BOT  generating  curve  defined  by  points  T^,  ...  with 

the  triangle  functions  centered  at  T^,  T ^ ,  etc.  Note  the  triangle  functions 
span  five  data  points,  with  adjacent  functions  overlapping  each  other. 


2.3  Input  Data  Description  and  Formats 

A  detailed  description  of  the  input  data  required  by  the  five  main-line 
MM/BOT  codes  is  provided  below.  READ  statements  and  formats  are  listed  in 
the  order  in  which  they  appear  within  the  programs,  followed  by  a  description 
of  the  required  data  for  each  of  the  main-line  MM/BOT  codes. 


This  data  set  Is  read  and  required  by  all  the  prograas. 

READ(5,1)NM0DE 
FORMAT (13) 

NMODE  -  Number  of  non-negative  modes  to  be  considered  (l.e. ,  there 
will  be  2*NM0DE-1  total  modes) . 

READ(5,2)NPT,NBAND 
FORMAT (213) 

tt  2t 

NPT  -  Number  of  diagonal  bands  to  be  used  In  each  Z_  .  Z  . 

U|ii 

Ztz  ,  and  Zzz  submatrix.  NPT*1  Indicates  that  only  the 
m,n  m,n 

diagonal  terms  are  nonzero  In  each  submatrix,  NPT*2  indi¬ 
cates  that  only  diagonal  and  off-diagonal  terms  are  non¬ 
zero,  etc.  If  NPT  >_  (NP-3)/2,  each  subnatrix  is  full, 
where  NP  is  described  below. 

NBAND  -  Number  of  submatrix  diagonal  bands  to  be  used  during  the 

inversion  of  ZgQ^.  NBAND* 1  indicates  that  only  the 

diagonal  Z  submatrices  are  to  be  used  during  inversion. 
m,n 

NBAND-2  indicates  that  diagonal  and  off-diagonal  Z^  n  sub¬ 
matrices  are  to  be  used,  etc.  If  NBAND  j>  2*NM0DE-1,  the 
entire  Z^  matrix  is  inverted. 

READ(5,3)NP,MC,BK 

FORMAT(2I3,E14.7) 

NP  -  Number  of  points  used  to  describe  the  BOT  generating  curve. 

(NP  must  be  odd.)  If  the  BOT  generating  curve  is  closed 
(l.e.,  the  first  and  last  points  coincide),  the  programs  will 
increase  NP  by  two  and  add  two  points  to  the  generating  curve 
(i.e.,  the  YH  and  XH  arrays  described  below).  This  new  NP 
should  be  used  in  all  definitions  involving  NP  (i.e.,  dimen¬ 
sions,  etc.). 

MC  -  Number  of  terms  to  be  used  in  the  numerical  integration  of  the 
Green's  function  [Equation  (34)  of  Volume  I]. 

BK  -  Wave  number  for  the  problem  (meters'1). 

READ(5,4)(YH(I),I-1,NP) 

F0RMAT(10F8.4) 

YH  -  Array  of  y  coordinates  for  the  generating  curve  (meters). 
READ(5,5)(XH(I),I-1,NP) 
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5  FORMAT (10F8. 4) 

XH  -  Array  of  x  coordinates  for  the  generating  curve  (neters). 
READ (5, 6) L 

6  FORMAT (F8. 4) 

L  -  Half  length  of  the  BOT  (meters;. 


2.3.2  Scattering  and  Radiation  Analysis  Data  Set 

This  data  set  is  read  by  the  BOTRA,  BOTSCM,  and  BOTSCB  programs.  Addi¬ 
tional  data  ate  required  by  some  of  the  above  programs  as  specified  later. 

READ(5,7)NANG,NT,PHII,THI 

7  FORMAT(2I3,2F8.4) 

NANG  -  Number  of  fixed  radiation  or  scattering  angles,  as  defined 

by  IPLANE  below.  NANG  radiation  or  scattering  patterns  will 
be  calculated. 

NT  -  Number  of  varied  radiation  or  scattering  angles,  as  defined 
by  IPLANE  below. 

PH11  -  <j>  angle  of  incident  wave  (degrees).  Used  only  in  BOTSCB 
program. 

THI  -  0  angle  of  incident  wave  (degrees).  Used  only  in  BOTSCB 
program. 

READ(5,8) (ANG(I) ,1-1 ,NANG) 

8  FORMAT (10F8. 4) 

ANG  -  Array  of  fixed  radiation  or  scattering  angles,  as  defined  by 
IPLANE  below. 

READ(5,9) (IPLANE(I) ,1-1, NANG) 

9  F0RMAT(10I8) 

IPLANE  -  Array  indicating  whether  the  corresponding  element  of 

array  ANG  is  a  $  or  6  angle.  IPLANE(I)-1  indicates  that 
ANG(I)  is  a  fixed  $  angle.  In  this  case,  $  is  fixed  at 
ANG(I)  and  0  varies  between  0°  and  180*  at  NT  equally 
spaced  angles.  (The  only  exception  is  program  BOTSCM, 
where  0  varies  between  0*  and  90s  at  NT  equally  spaced 
angles.)  IPLANE (I)-2  indicates  that  ANG(I)  is  a  fixed  0 
angle.  In  this  case,  6  is  fixed  at  ANG(I)  and  $  varies 
between  0°  and  180°. 


2.3.3  Antenna  Related  Data  Set 


This  data  set  la  read  only  by  the  BOTRA  program  and  should  be  deleted 
when  running  the  BOTSCB  program. 

BEAD(5a10)NSA 

10  FORMAT  (13) 

NSA  -  Number  of  slot  antennas  on  the  BOT. 

READ(5,11)(IS(K) ,K-1,NSA) 

11  FORMAT(10I8) 

IS (K)  -  Triangle  fiuctlon  at  which  slot  antenna  K  is  located 
(centered) . 

READ(5,12) (Z0(K) ,K-1,NSA) 

12  FORMAT (10F8. 4) 

Z0(K)  -  Starting  Z  coordinate  for  antenna  K  (meters). 

READ(5,I3) (Z1(K) ,K-1,NSA) 

13  FORMAT(10F8.4) 

Z1(K)  -  Ending  Z  coordinate  for  antenna  K  (meters).  (ZO(K)  < 

Z1(K)  for  all  K.) 

READ(5,14) (E0(K) ,K-1,NSA) 

14  FORMAT (10F8. 4) 

E0(K)  -  Constant  excitation  across  slot  antenna  K  [Equation  (59) 
of  Volume  I].  (E0(K)  is  complex.) 

READ(5,15) (TEXC(K) ,K-1,NSA) 

15  FORMAT (10F8. 4) 

TEXC(K)  -  Indicates  t  excitation  on  slot  antenna  K  when  TEXC(K)  i  0. 
This  is  the  (U£)^  term  in  Equation  (47)  of  Volume  I. 

READ (5, 16) (ZEXC(K) ,K-1,NSA) 

16  FORMAT (10F8. 4) 

ZEXC(K)  -  Indicates  *  excitation  on  slot  antenna  K  when  ZEXC(K)  +  0. 
This  is  the  (U*)^  term  in  Equation  (47)  of  Volume  I. 

2.3.4  Aperture  Analysis  Data  Set 

This  data  set  is  read  only  by  the  BOTSCB  program  and  should  be  deleted 
when  running  the  BOTRA  program. 

READ(5,17)Z0,Z1,Y0,X0,Y1,X1 


■  {if? 

■  "X  ■ 

5:  l 

,  .w 


17  F0RMAT(6F8.4) 

ZO  -  Starting  z  coordinate  for  the  aperture  (meters). 

Z1  -  Ending  z  coordinate  for  the  aperture  (meters).  If  Z1  £  ZO,  an 
aperture  is  not  present,  and  the  remaining  parameters  in  this 
data  set  are  ignored. 

YO  -  Starting  y  coordinate  for  the  aperture  (meters) . 

XO  -  Starting  x  coordinate  for  the  aperture  (meters).  If  the  start 

of  the  aperture  is  on  generating  curve  segment  I,  the  point 
(XO,YO)  should  be  approximately  on  the  line  segment  joining 
(XH(I).YH(I))  and  (XH(I+1) ,YH(I+1)) .  To  be  specific,  XO  must 
lie  between  XH(I)  and  XH(I+1) ,  and  YO  must  lie  between  YH(I) 
and  YH(I+1). 

Y1  -  Ending  y  coordinate  for  the  aperture  (meters) . 

XI  -  Ending  x  coordinate  for  the  aperture  (meters).  If  the  end  of 

the  aperture  is  on  generating  curve  segment  J,  the  point 
(XI, Yl)  must  satisfy  the  same  conditions  as  (XO,YO)  above 
with  J  replacing  I.  The  aperture  is  assumed  to  start  at 
(XO,YO)  and  extend  to  (XI, Yl)  in  the  t  direction  (i.e.,  in 
the  direction  of  the  generating  curve).  If  the  aperture  does 
not  cover  the  peaks  of  at  least  two  triangle  functions,  the 
program  will  print  an  error  message  and  stop  execution. 

2.3.5  Near-Field  Analysis  Data  Set 

This  data  set  is  read  by  the  BOTRA  and  BOTSCB  programs. 

READ(5,18)NTEST 

18  FORMAT (13) 

NTEST  -  Number  of  test  points  at  which  near-field  radiation  or 

scattering  is  to  be  calculated.  NTEST  may  be  set  to  zero. 
Repeat  the  following  NTEST  times: 

READ(5 , 19) ZTEST , YTEST , XTEST 

19  FORMAT(3F8.4) 

ZTEST  -  z  coordinate  of  the  test  point  (meters). 

YTEST  -  y  coordinate  of  the  test  point  (meters). 

XTEST  -  x  coordinate  of  the  test  point  (meters) . 
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2.4  Parameter  Selection 

The  choice  of  moat  of  the  parameters  In  Section  2.3  is  specified  by  the 
user  depending  upon  the  explicit  requirements  of  the  problem.  For  example, 
NSA  is  set  by  the  number  of  slot  antennas  on  the  BOT.  Similarly,  BK  is 
determined  by  the  frequency  at  which  the  MM/ BOT  analysis  is  carried  out.  On 
the  other  hand,  the  choice  of  the  parameters  NP,  MC,  and  NMODE  is  based  upon 
the  requirements  of  the  MM/BOT  theory,  as  explained  below. 

The  parameter  NP  specifies  the  number  of  coordinate  points  used  to 
describe  the  BOT  generating  curve  in  the  t-direction.  These  points  in  turn 
define  the  segmentation  (strips)  of  the  BOT  surface.  In  general,  the  sepa¬ 
ration  of  the  points  should  not  exceed  0.15  X.  If  the  BOT  has  a  rapidly 
changing  surface  geometry  with  sharp  curvatures,  a  higher  density  of  points 
may  be  necessary.  The  resulting  number  of  triangle  functions  used  to  expand 
the  surface  currents  on  the  BOT  is  determined  by  NM  *  (NP-3)/2.  Since  at 
least  one  half  of  a  triangle  function  should  subtend  an  active  aperture 
antenna  in  BOTRA  and  at  least  two  triangle  functions  must  intercept  a  passive 
aperture  in  BOTSCB,  the  value  of  NP  is  specified  by  these  conditions  as  well 
as  the  strip  width.  As  a  general  principle,  the  accuracy  of  the  analysis  is 
improved  by  increasing  NP. 

The  parameter  MC  is  used  in  the  numerical  evaluation  of  the  integrated 
Green's  kernel  [Equation  (26)  in  Volume  I).  A  necessary  condition  for 
the  approximations  used  in  obtaining  G  is  that  MC  >  8(NM0DE-1) . 

The  parameter  NMODE  is  set  by  the  length  of  the  BOT.  In  general,  the 
minim, im  requirement  is  that  NMODE  >  2L/X,  where  L  is  the  (axial)  half-length 
of  the  BOT.  This  requirement  is  comparable  to  the  MM/BOR  analysis  require¬ 
ment  that  the  maximum  circumferential  modes  used  be  n  ^  ttD/X,  where  D  is  the 
largest  diameter  of  the  BOR.  As  a  general  observation,  the  accuracy  of  the 
analysis  increases  and  the  spatial  resolution  of  the  surface  currents  on  the 
BOT  is  improved  as  NMODE  is  increased.  This  is  particularly  true  for  the 
edge  currents.  However,  practical  computer  main  memory  limitations  usually 
set  the  upper  limit  on  NMODE. 


2.5  Program  Dimensions 


The  MM/BOT  programs  (Appendix  C)  are  currently  set  up  to  handle  problems 
with  the  following  set  of  parameters.  (Some  of  the  arrays  are  overdimen¬ 
sioned  In  the  listings.) 

NMODE  <  4 
NBAND  <  2*NM0DE-1 

MP  <_  19  (17  for  a  closed  generating  curve) 

MC  <  50 
NANG  <  6, 

NT  <  91 
NSA  <  20 

NAM  12  (NAM  is  the  number  of  triangle  function  peaks  subtended  by  an 
aperture) . 

For  a  different  set  of  input  parameters,  the  minimum  dimensions  required 
for  each  program  are  listed  below.  Table  1  contains  definitions  of  the 
parameters  used  as  the  subscripts  in  the  arrays  enumerated  in  the  dimension 
statements. 


TABLE  1.  DEFINITION  OF  DIMENSION  STATEMENT  INDICES 


Parameter 

Definition 

LS 

NP-3 

MC 

Input 

NAM 

Number  of  triangle  function  peaks  subtended  by 
an  aperture 

NANG 

Input 

NM 

(NP-31/2 

NMODE 

Input 

NP 

Number  of  points  on  the  generating  curve  (input). 

If  the  curve  is  closed,  use  NP+2  in  place  of  NP  in 
all  definitions. 

NSA 

Input 

NT 

Input 

OPTf-OMI-lOf 


BOTZSS  Minimum  Dimensions 

COMPLEX  Z(LS*LS) ,G((NP-l)*NP/2) 

DIMENSION  YH(NP) ,XH(NP) ,DH(NP-1) 

DIMENSION  SV(NP-l) ,CV(NP-1) ,XS(NP-1) ,YS(NP-1) 


DIMENSION  UMN(MC) , JK(4) 

DIMENSION  TP(4*NM) ,T(4*NM)  ,TZ(4*NM) 
BOTINV  Minimum  Dlnenalons 

DIMENSION  YH(NP) ,XH(NP) 

COMPLEX  Z(K1) ,ZI(K2) ,WORK(K3) 
DIMENSION  NZ (K4) 

COMMON  NM,JK(4),LR(K5) 


where  K1 

through  K5  depend  on  NBAND  and  are  defined 

below. 

NBAND 

1 

<2*NM0DE-1 

>2*NMODE-l 

K1 

0 

LS2*  { (2*NM0DE— 1) 

* (2*NBAND-1) -(NBAND- 1) *NBAND } 

|LS*(2*NM0DE-1) } 

K2 

LS2 

LS2* ( 2*NM0DE-1) 

LS2 

K3 

0 

LS 

0 

K4 

0 

2*NM0DE-1 

0 

K5 

LS 

LS 

LS* (2*NMODE-l) 

BOTRA  Minimum  Dimensions 


DIMENSION  IS(NSA) ,ZO(NSA) ,Z1(NSA) .TEXC(NSA) ,ZEXC(NSA) 

COMPLEX  EO (NS A) 

COMPLEX  VM(LS) ,VN(LS) 

COMPLEX  GT(NT) ,GP(NT) 

COMPLEX  ESC (3) ,HSC(3) ,CUR(LS*(2*NM0DE-1)) 

COMPLEX  Y(LS*LS) ,RT(LS) ,RP(LS) 

COMMON / BODY / DH (NP-1 ) ,R(NP-1) ,SP(NP-1) ,CP(NP-1) ,SV(NP-1) ,CV(NP-1) , 
BK,L,NP,T(4*NM) ,TZ(4*NM) 

COMMON/TEST /MC,  GAMMA,  TP  (4*NM)  ,XS(NP-1)  ,YS(NP-1) 

DIMENSION  YH(NP) ,]ffl(NP) ,THR(NT) ,PHIR(NT) 

DIMENSION  ANG(NANG) , IPLANE (NANG) 

Subroutine  PLANE  dimensions  are  as  follows: 

COMMON/BODY / DH (NP-1 ) ,R(NP-1) ,SP(NP-1) ,CP(NP-1) ,SV(NP-1) ,CV(NP-1) , 
BK,L,NP,T(4*NM) ,TZ(4*NM) 

Subroutine  NEARB  dimensions  are  as  follows: 

COMPLEX  G(NP-l) ,HO(NP-l) ,H1(NP-1) 

COMMON/BODY/DH(NP-1) ,R(NP-1) ,SP(NP-1) ,CP(NP-1) ,SV(NP-1) ,CV(NP-1) , 
BK,L,NP,T(4*NM) ,TZ(4*NM) 

COMMON / TEST / MC , GAMMA , TP  ( 4*NM)  ,XS(NP-1)  ,YS(NP-1) 
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BOTSCM  Minimum  Dimensions 

COMPLEX  STT(NT) ,SPP(NT) ,STP(NT) ,SPT(NT) 

COMPLEX  YT(LS) ,YP(LS) 

COMPLEX  Y(LS*LS),RT(NT*LS),RP(NT*LS) 

COMMON /BODY / DH (NP-1) ,R(NP-1) ,SP(NP-1) ,CP(NP-1) ,SV(NP-1) ,CV(NP-1) , 
BK,L,NP,T(4*NM) ,TZ(4*NM) 

DIMENSION  YH(NP) ,XH(NP) ,THS(NT) ,PHIS(NT) 

DIMENSION  ANG(NANG) .IPLANE(NANG) 

Subroutine  PLANE  dimensions  are  given  under  BOTRA  dimensions. 

BOTSCB  Minimum  Dimensions 

COMPLEX  STT (NT) , SPP (NT) , STP (NT) , SPT (NT) 

COMPLEX  ESCT(3) ,HSCT(3) ,CURT(LS*(2*NMODE-l); 

COMPLEX  ESCP(3) ,HSCP(3) ,CURP(LS*(2*NM0DE-1)) 

COMPLEX  Y(LS*LS) ,RT(LS) ,RP(LS) 

COMMON/BODY/DH(NP-1) ,R(NP-1) ,SP (NP-1) , CP (NP-1) ,SV(NP-1) ,CV(NP-1) , 
BK,L,NP,T(4*NM) ,TZ(4*NM) 

COMMON/TEST/MC, GAMMA, TP(4*NM) ,XS (NP-1) ,YS(NP-1) 

COMMON/S/LT  (NAM)  ,VLOW(NAM)  .VHGH(NAM)  ,  IT  (NAM) 

DIMENSION  YH(NP) ,XH(NP) ,THS(NT) ,PHIS(NT) 

DIMENSION  ANG(NANG),IPLANE (NANG) 

Subroutine  PLANE  dimensions  are  given  under  BOTRA  dimensions. 
Subroutine  NEARB  dimensions  are  given  under  BOTRA  dimensions. 
Subroutine  ASYMAC  dimensions  are  as  follows: 

COMPLEX  YA(NAM*NAM) ,CAT(2*NAM) ,CAP(2*NAM) ,EVT(2*NAM) ,EVP(2*NAM) 
DIMENSION  FPP (3 ,NAM) FPPZ (3 ,NAM) ,LR(2*NAM) 

COMMON/S/LT (NAM) ,VLOW(NAM) .VHGH(NAM) , IT (NAM) 

Subroutine  APPAR  dimensions  are  as  follows: 

COMMON/S/LT (NAM) .VLOW(NAM) , VHGH (NAM) , IT (NAM) 


2.6  Sample  Problems 

.In  this  secCion,  four  sample  problems  are  considered  to  illustrate  the 
use  of  the  MM/BOT  algorithm.  Sample  problem  1  exercises  all  the  main-line 
programs  given  in  Figure  1.  The  inputs,  outputs,  and  selected  intermediate 
results  are  reproduced  here  to  provide  a  check  case  for  the  proper  function- 
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ing  of  Che  codes.  Problem  2  demonstrates  the  radiation  analysis  for  aperture 
antennas  embedded  in  an  asymmetric  BOT  (i.e.,  wing  section).  Problem  3  con¬ 
siders  the  scattering  from  a  square  cylinder  of  finite  length.  Finally, 
problem  4  demonstrates  the  use  of  the  codes  to  compute  the  aperture-coupled 
electric  and  magnetic  fields  inside  a  BOX,  approximating  a  fuselage  equipment 
compartment  when  subject  to  electromagnetic  illumination. 

2.6.1  Problem  la 

Consider  a  right-circular  cylinder  of  2.76  X  length  and  0.216  X  radius, 
with  an  embedded  ^-polarized  aperture  antenna  at  <{>  *  0°.  The  aperture  is 
fed  uniformly,  subtends  a  45°  opening,  and  is  2.06  X  long  in  the  axial  direc¬ 
tion  (see  Figure  3a). 

a)  Calculate  the  power  gain  patterns  in  the  horizontal  ($  *  0,  180°) 
plane  and  the  roll  (8  ■  +90°)  planes  in  the  6  and  <p  polarizations. 

b)  Compute  the  currents  on  the  cylinder  surface. 

c)  Compute  the  electric  and  magnetic  field  components  (near  fields)  at 
test  points  on  a  line  bisecting  the  aperture,  i.e.,  at  XTEST  ■ 
14.256,  19.44,  71.28,  150,  and  300  m  where  the  corresponding 
YTEST  -  ZTEST  -  0.0. 


Figure  3a.  Slotted  cylinder  for  problem  1 . 


Solution  -  The  problem  was  solved  by  running  the  programs  BOTSEG,  BOTZSS, 
BOTINV,  and  BOTRA  in  sequence.  The  calculations  were  carried  out  at  10  MHz 
(X  ■  30  m).  The  cylinder  was  represented  by  NP  ■  17  points  around  the  cir¬ 
cumference. 

Execution  of  BOTSEG  is  shown  in  Figure  3b,  together  with  the  data  file 
for  the  coordinates  of  the  problem  and  a  plot  of  the  HP  data  points  and 
triangle  functions  for  the  cylinder  cross  section.  (Note  Figure  3b  was 
generated  in  a,  time-share  mode.  The  question  marks  at  the  beginning  of  each 
line  are  prompt  signs  used  in  the  time-share  mode.)  The  generated  file  was 
used  as  the  basic  input  file  for  the  programs  BOTZSS,  BOTINV,  and  BOTRA.  Four 
modes  were  used  in  the  calculations,  i.e.,  NMODE  *  4.  Partial  listings  of  the 
inputs  and  outputs  from  these  programs  are  shown  in  the  subsequent  figures. 
Specifically,  the  input  data  for  execution  of  the  programs  are  given  in  Figure 
4.  For  reference,  the  variable  of  the  data  set  are  labeled  with  the  appli¬ 
cable  format  statements  in  parenthesis.  Note  the  aperture  coincided  with  the 
fourth  triangle  function  (IS  *  4)  and  only  one  triangle  function  was  used  to 
span  the  aperture  (NSA  -1).  If  a  non-uniform  aperture  excitation  is  desired, 
then  more  triangle  functions  should  be  used  to  span  the  aperture  each  with  a 
different  Eq.  Partial  outputs  from  BOTZSS,  BOTINV,  and  BOTRA  are  shown  in 
Figures  5,  6,  and  7,  respectively.  The  radiation  power  gain  for  the  slotted 
cylinder  for  the  horizontal  and  roll  planes  normalized  to  an  isotropic  radia¬ 
tor  is  summarized  in  Figure  7b.  (The  comparison  of  these  results  with  the 
MM/BOR  analysis  is  given  in  Figure  10  of  Volume  I.)  Partial  output  of  the 
currents  on  the  cylinder  is  plotted  in  Figure  7c.  The  electric  and  magnetic 
fields  computed  in  the  near-field  analysis  are  listed  in  Figure  7d.  (For  fur¬ 
ther  discussion  of  these  results,  see  Section  8  of  Volume  I.) 

The  foregoing  calculations  were  carried  out  at  10  MHz.  If  the  dimen¬ 
sions  of  the  body  (BOT)  and  the  antenna  are  given  initially  in  terms  of 
wavelength,  any  convenient  frequency  can  be  chosen  in  the  setup  procedure  for 
carrying  out  the  computations.  If  the  data  are  to  be  compared  with  range 
measurements  at  a  given  frequency,  for  ease  of  data  interpretation,  the 
calculations  are  also  done  at  that  frequency. 
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7  0.0  0.0 
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Figure  3b.  Execution  of  BOTSEG  for  problem  1 ,  with  triengle  functions  numbered  on  the  plot. 
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Figure  5a.  Partial  output  of  BOTZSS  for  problem  la. 


L1&  V:**-**  .if- 


«  H  H  •  O  •  H 

s  s  s  1 1 1  a 


M  «  N  « 

r  t  t  r 

m  S„i**5*kw  »ms«» 


U  :  s  s  i  i  s 


•  «4  # 

:u 


•  •ill 


?  •  i  7  m  ill# 


•  «««S»9«*I  *4 


?  ? 

m  m 

9  7  9  S  * 


«  •  M  • 

4  »  «  7 


T  •  T  T 

•  O  m  *4  a*  ♦  •4«*S««*«|£  «  4  H  » 

Si  #  •  m  •  ><  «  «t  e 

•  ►  (K  »  K  e«««4AOO  0««t4«40i  «  «  »  • 

#•«!••  •  ••■»•••■*  •  •  •  • 
«*  «•  *1  A*  *>  •>  M  M  M  o  M  A*  M  •  ******  *******  •  •  4*  <  41 

III  lilt  III  I  I  I  I  I 

«  M  i  2  •*  *» 

T*4  O  O  0*4 

•  T  •  it 

Ml  |g  in  Mi  Ml  III 

«  A  «  #  >4  4  »«OB04»4  MAN* 

«•«»•«  A4AAA4M4  MN»N»NM«  O  H  4  4 

A  •  A  •  M  O  A4AOA4MO  4AAOAM44  4  4  4  4 

**  *  N  -  •  m  •  t  t  f  r  *  i  •  i*  •  i  N 


•4  «  <•  At 

A  4  A  N 


?  7  7  ?  T  7 


stSSSsssssS?  sssssssgssssssss  ssssssssssssssss  swsssss  ssssss 

•  •  •  40  •  •«••••••  *0  •  •  •  *M  ••»«•••  *A  •  •  ••  ««4«  •*!  •  •  •  «4*  •  «A  ••••«•!•  •••••••• 

-rrrr—^-rr  7-7- r^T?  mr:  r-r— rr-r-Tr*  ttttttt?  rrrruT? 

4  A  A  4  4*4  A  A4  A  A4 

t  t  r  ?  2iS:!! 

44444A444AA444AA  A40AAA0AA  44AA4  4»  400 4><4  Of  44044004  A|i4>4WM| 

44AA4 mmmmmmm  »4i4N4i  4«»i  aop  *i»»»i«N»ONN»»4ia  — 114444044  1  •o^4f«4 

4AA44  4AA0^<4l44i04W4»  4AH4SAA»4^<4g<4  *0  AAAA444WA44W44  ••  44444444  4444A4  44 

•  •••••••••••••••  •••  *4A  *A  »4h  m4  »■  4AI  A  •  «AAA  >AAAA»  4  «AH  ••••••44  ••#••••• 

P4WNI 114  441  «■  •  44>I4A4A  HAAA  •  «A  M  •  *4  *A  •  •  »A4  •  •  *4  •  •  •  •  »  M  •  •  A4»4>A4W  I  II*  »W»  *4 


?  t  r?  ?  r  r  ? 

Ml  Mi  U  HI  Uf  Ml  Ml  »*# 

4aom«4a>4<i  t+m'i&tm  m*  ■  *>i+m*+*t+*  9mt*5**+  44444441  n  iBWuft;  **++>+**+++  m*mtm*t** 

areggsggsatgnsg  ssssKssssssssss  s^tsittt^ststs  sasjsss  sxsssss 

•  0O  •••••••••  «0  •  •  •  *4  •  (A  ••••4«  »N  •  •  •  4A  •  «4  •  •  •  *4  •  *4  •  •  M  •  •  •  •  •  •  •  •  •*  •  •  •  • 

«4l4UAi  W»A»I  «4»H*  i<AN  NA  ••Al  «4«»44At  *NN  «AA  44  *44*4444  *44  144  41  44  f4l>4  g 

III  I  I  •  I  I  I  I  I  MM  »  •  I  •  M  I  IITTIMI  TTTiTTm 

$  ~  sj  88  S  83  2  33 

'  '  'i  AAiiliii 

SXSSSSXSSSSSSX99  wtgfffa  85848888 

MOWA04MW40AW04A  44Q44W44444>Mf44  44444|4444AM444  44444444  4  *A»»»  W 

#••••#••••••••••  •  *44  *4  «44  4«  <44  •  •  «A4P  *044444  MN4  •  ••••••••  •••••••• 

li  W4#<4#i4»##< •—+14*441  AM 4  •  «A  M  •  «4  *A  •  •*•«•  4  *  •  »4  •  •  •  •  •  M  •  »  *4  44444444  ■•  IW44|4» 

•  Mil  I  ••  •  •  •  I  •  ••  III  III 

4  4  N  4  4  4  *•*» 

?  r  ?  ?  rr  rt 

*£9t£SSSX9mSS89m  y*8^S**5^8#?85SS  S2233222  *£«$£££$ 


•  •••••••••*•••••  A  •  *4  •  •  •  •  *4  •  •  •  •  *4  IjA  •  •  •  *N  ••#»••••  ,M  ««.**•.*  »««.».»? 

rr*"rrr  r-.-rTTP.-rr*  —  rrr’rrTT*  tttttttt  Trrrrrrr 


il  n  siaxg 

*****************  ********** ******  jMggMgWg 


JS'W.SttiTCSwW  88*iwKS^Sw  ¥?£&S  W 
r«««N«uwM  « -rr *r •  — r*"r  -mVt.vv-vtv  *~7~r~rr  r? 


Partial  output  of  BOTZSS  -  impadanoa  matrices  for  problem 
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Figura  Ba.  Partial  output  of  BOTINV  for  probfem  1*. 
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Admittance  matrices  for  mode  numbers  M,  N  are  output  by 
columns.  In  output  below,  Y1  in  a  complex  8x8  matrix. 


Partial  output  of  BOTINV  •  admittance  matrica*  for  m  *  n  *  -3  for  problem  la. 


Fifin  7a.  Partial  output  of  BOTR A  for  problem  la. 


Figure  7a.  Partial  output  of  BOTRA  for  problem  la.  (concluded) 
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Figure  7b.  Partial  output  of  BOTRA  for  proMam  la. 
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Figure  7b.  Partial  output  of  BOTRA  for  problem  ta.  (continued) 
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Figure  7b.  Partial  output  of  BOTRA  for  problam  la.  (continuad) 
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Figure  7b.  Partial  output  of  BOTH  A  for  problam  la.  (oondudad) 
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Figure  7c.  Partial  output  of  BOTRA  for  probiam  la. 
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Figure  7c.  Mai  output  of  0OTRA  for  problem  la.  (concluded) 


2.6.2  Problem  lb 


Compute  the  monostatic  (baclcscatter)  cross  sections  at  10  Mis  for  the 
following  received  and  transmitted  polarizations  06,  4><f>,  <f>8,  and  8$  in  the 
horizontal  ($  *  0)  and  roll  (6  *  90)  planes  for  the  cylinder  configuration 
used  in  problem  la.  For  this  part  of  the  problem,  assume  that  there  is  no 
aperture  on  the  cylinder.  Compute  the  cross  sections  at  46 (NT)  equi- 
spaced  angles. 

Solution  -  Since  the  basic  body  configuration  was  not  altered,  the  results 
from  BOTZSS  and  hence  BOTINV  in  problem  la  are  used.  The  cross  sections  were 
obtained  from  BOTSCM.  Figure  8  depicts  the  input  data  to  execute  BOTSCM, 
Partial  listings  of  the  output  from  BOTSCM  including  the  scattering  cross 
sections  are  given  in  Figure  9. 
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Figure  8.  Input  data  for  execution  of  BOTSCM  for  problem  1b. 
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Figura  9a.  Partial  output  of  BOTSCM  for  probiam  1b. 


2.6.3  Problem  lc 

Consider  the  cylinder  in  problem  la  with  a  passive  aperture  subtending 
an  angle  of  45°  and  extending  from  -27.9  <_  Z  27.9  m  along  the  length  of  the 
BOT.  Let  the  edges  of  the  aperture  parallel  to  the  z-axis  be  defined  by  the 
coordinates  YO  -  4.5822,  XO  -  11.062,  and  Y1  -  -0.001,  XI  -  12.95.  (Let  NP  - 
17  and  NMODE  -  4.) 

a)  Calculate  the  bistatic  cross  sections  when  a  TM  (0-polarized)  wave 
illuminates  the  body  broadside  (8  »  90°) .  Repeat  the  analysis  for  a 
TE  ((((-polarized)  wave.  For  this  calculation,  assume  that  the  body 
does  not  have  an  aperture. 

b)  Calculate  the  electric  and  magnetic  fields  penetrating  the  aperture 

when  the  body  is  again  illuminated  broadside  with  a  10  MHz  TM  or  TE 
wave,  along  a  line  bisecting  the  aperture  (i.e.,  the  incident  angles 
for  the  wave  are  0^  *  90°,  ■  22.5°).  Sample  the  fields  at  the 

following  points,  measured  in  meters,  within  the  cylinder  along  a 
line  bisecting  the  aperture: 
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YTEST 

XTEST 
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12.4667 

i  following  points  outside 

the  aperture: 

0.0 

0.0 

14.2560 

0.0 

0.0 

19.44 

0.0 

0.0 

71.28 

0.0 

0.0 

150 

0.0 

0.0 

300 

Solution  -  The  solution  to  the  two  parts  of  problem  lc  was  obtained  by 
running  BOTSCB,  using  BOTINV  output  from  problem  la  as  the  input.  In  this 
analysis,  the  aperture  was  intercepted  by  two  triangle  functions.  (For 
larger  aperture  angles,  more  triangle  functions  fall  within  the  aperture 
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opening.)  Figures  10  and  11  depict  the  coordinate  file  and  partial  listings 
of  the  inputs  and  outputs  of  BOTSCB.  The  bistatic  scattering  cross  section 
for  the  apertureless  body  and  the  corresponding  currents  (partial  listing) 
for  6  and  4>-polarized  incident  waves  are  given  in  Figures  lib  and  11c, 
respectively.  These  outputs  constitute  the  solution  of  the  first  part  of 
problem  lc.  The  solution  of  the  second  part  of  the  problem,  involving  the 
cylinder  with  the  passive  aperture,  is  summarized  in  Figure  lid. 
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Figure  10.  Input  data  for  execution  of  BOTSCB  for  problem  1c. 


39 


,  f  ft*  .■•' 

ilinriirfia  m  ii  r(n- 


NFT 

N8AN0 

NNODE 

NF 

NC 

20 

14 

4 

17 

50 

TH 

-im 

-ItoBOT  — *IS©  -3 19867 

6*4800 

-4.5821 

XH 

xi-m 

unit 

iMi 

v.m 

HALF-LENGTH  of  «OT  •  41.4000 

NUMBER  OF  FIXED  ANGLES  •  2 

NUMBER  OF  ANGLES  PER  FIXED  ANGLE  ■  46 


BK 

.  2094 39 7 E 400 

2-*r,#  °*oood  -1*4™ 
®*251S  1is||||  12**667  12*9600  12*4667 


FIXED  ANGLE 


22.5 


90.0 


CODE  I  PHI  FIXED  ■  1*  THETA  FIXED  *2) 

1 

2 


N*  -3  N«  -3 
T 

•9070311E-04  — .47D10B7E-03  .6824476E-04  -. 5931252F-C3 

H«  -2  M*  -3 
T 

-• 7578504E— 04  .1106460E-04  -.5825230E-04  .2733491E-04 

N*  -1  N*  -3 
T 

.88183036-04  . 12 7903BE-04  .62791096-04  -.26567026-04 

N*  0  N*  -3 
T 

-•1079845E-03  -.3263554E-04  -.69704426-04  .28628486-04 


Fiaura  11a.  Partial  output  of  B0TSC8  for  proMam  1c. 
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Figure  11b.  Partial  output  of  BOTSCB  for  problem  1c.  (continued) 


K  -  •*  '! 


44 


g  8 

t  ’=2 
□  IS  <0  .£ 
o  C  o  ■*- 
■P.S4- 

N  I  “-e-J 

■o  G  £  •®-S 

5  |  8  5  I 

^  *D  o  *n- 


M  Am  AA  ##  NN  MM  IMA 

TTTTTTff??TfffffTT?? 


II  II  II 


AO  AA  Ml  AA  Ml  A  A  Ml  MM  Am  #•*■  Ml  AA  AA  AA 

oo  oo  oo  oo  oo  OO  OO  o©  oo  oo  oo  oo  op 

TT  II  II  II  II  TT  i  •  •  •  I  i  TT  XT  II  Tl  •• 


AO  AM 

AO  AM 

••  AM 

AO  AA 


«*  M  AO  OO  H»  AO  AO  ON  CN  AO  AM 

»«  AM  OO  MA  AA  AM  oo  »A  OO  AA  AM 

N-A  AA  MO  M  M  AA  HO  AO  N-A  AA  AM 

nv  MM.  MA  AM  AA  **M  M  MM  AA  AA  AA 


mm  AA  MM  MM 

OO  OO  OO  MM 

•  I  II  II  II 


mm  mm  OM  Am 


AA  AA  Am  AA 

oo  oo  oo  oo 

•  I  II  II  II 


mm  AA  MA  MM 

op  oo  oo  mm 

TT  ii  ii  ii 


AO  AA 

•  •  •  • 

I  I 


MM  MA  AN- 


1 

• 

1 

1 

6 

A* 

E 

93 

z 

AA 

mm 

AA 

mm 

AA 

o 

OO 

1  1 

OO 

1 « 

7? 

?? 

oe 

1  • 

a 

A. 

OO 

AA 

UNM 

AA 

OA 

MA 

cm 

MM 

UNM 

Mm 

OA 

MO 

mM 

UNM 

AM 

AM 

MA 

w 

& 

CD 

•  • 

•  • 

•  • 

•  • 

•  • 

O 

1 

• 

N 

AA 

mm 

AA 

AA 

AA 

O 

OO 

OO 

OO 

OO 

OO 

CD 

UN 

UIUI 

ss 

MO 

ON. 
OA 
•  « 

mo 

OO 

NA 

mo 

A® 

mm 

mM 

AM 

AM 

O 

3 

MA 

MM 

mo 

mo 

AA 

a. 

1 

• 

1 

1 

1 

M 

3 

O 

AA 

mm 

mA 

OO 

oo 

oo 

•  1 

•  l 

1  l 

Nui 

U.U 

UIUI 

AA 

m© 

OO 

OA 

AO 

AO 

AN. 

ON- 

Mm 

MA 

MM 

mM 

• 

H 

AA 

AA 

AA 

OO 

1  1 

?? 

OO 

1  1 

UNM 

UiUl 

UiU. 

MO 

oo 

OA 

MA 

OM 

ON- 

AM 

AM 

MA 

OO 

MM 

m  • 

M  |l 

AA  A  A  MM  MM  AA  A  A  AA  AA  TlA  AA 
OO  OO  OO  OO  OO  OO  OO  OO  OO 

II  II  II  II  II  II  II  II  II  II 

UM  Ultil  UU'  UJUi  tiW  U.UI  UiUl  IM'  UM)  Ibli) 

AA  AA  MA  AA  AA  AA  OO  iXD  AA  MO 

AN-  AN.  MO  AO  AN-  AM  OO  MA  AA  OO 

MM  AA  MO  MA  MM  AA  MA  AM  IN  ON* 

AM  AA  AM  MM  AM  AA  AM  M9  MA  NN 


AA  AA 

??  ?? 


-A  AO  ®M  MA  MO  AO  OO  AM  MO 

MA  MO  MN  AO  AA  MO  OM  MA  OA 

MO  AA  MA  NC  MO  AM  MO  ON 
M  AO  Mm  04  mA  AO  MM  MA 


I  O  I  M  I  M  I 


cop  Omm 
u.  I  T  ii  i  i 

WU«  UN' 

A  AM  A  AO 
MOM  MO  A 
OOM  *  AM 
lifNA  Ui  MM 


coo  ope  coo  ooo  ooo 

u.  •  i  u. TT  u.  i  i  u. Tl  u.  i  i 

imj  mu  umr  UNN 

AAA  A  AO  A  AO  ACC  A  AC*M 


oop  ooo  opo 

IV  I  |  Ui  I  I  lii  T  | 

NliN.  M  UMi  MuiUi 
OM#  UM<  O  OA 
IN#M  UiOA  UlOM 
• OM  KBN  N OM 
AMU  MAM  MOA 
O  •  •  O  •  •  On* 
•  III  II 


•  II  «  I 

3  9 

o  o 

MM  MM 
UUO  OMm 
Ui  |  I  u>  •  I 

MIMA'  NUW 
UHV  O AM 
MINM  UIOA 
AAA  A AM 


MOA  MMM 

o  •  •  o  •  • 
I  I  III 


AA  AA  AA 

coo  coo  ooo 

Mill  Ui  I  I  Ui  |  I 


Ui  AN-  Ui  AO  UI  AM 
■  OM  AMO  AMA 
mmm  MAM  MOA 
O  •  •  0*1  pi* 

ii  i  T  I 


Plot  of  currents  along  BOT  for 
^-polarized  incident  field 


T  DIRECTED  CURRENTS  ON  TRIANGLE  FUNCTION  1 
♦  MAGNITUDE 


•  135  £-02 
•134E-02 
•133E-02 

•  133E-02 
•132E-Q2 
•131§-02 
•130E-02 
•129E-02 
.1286-02 
•128E-02 
•127E-02 
•126E-02 
•125E-02 
•12AE-02 
•123E-02 
•123E-02 
• 1226-02 
.1216-02 
•120E-02 
.119E-02 


e« 


•  « 


*  * 

*  *  * 


.1196-02  I 


*  e 
* 


*  ♦ 
* 


♦  PHASE 

180*0 

162.0 

1AA.0 

126.0 

106.0 

90.0 

ih 

o.o 

-18.0 

:k:8 

-II:S 

:l«:8 

-14A.0 

-162.0 

-180.0 


I 

I 

■~i  *  17“ 

-1.0 

-.50 

Z/L 

0.  .50 

1. 

Z-DIRECTED  CURRENTS 

ON  TRIANGLE 

FUNCTION  1 

•  MAGNITUDE 

♦  PHASE 


.321E-03 

•305E-03 

.2996-03 

.272E-03 

•256F-03 

.2506-03 

.225F-03 

•209E-03 

•160E-03 

.1556-03 

.1286-03 

•801E-05 

• 139E-11 


*****!*&&&*'.*#*■ 


a?ertu*c  analysis 


•27?$000  27?&00C  4l$B22  11?$620  -!oC10  lzfJj^O  Aperture  coordinate 


LT  VL  Ou 


.9001 

2.9204 


2.9204  -1 

.0014  1 


IT  IT 


4  ~i  alioieco  llrsif eSS  :HMi:  t?  rlr^ieoc  :mv» 

“hsflttl  ziiinm  ::$SKf4i  Y?  ' 

'littfclt  tftfcit  ziittmt  Yart 

»»s  r Aperture 

'•f 

T44 

-:i2!8§:8?  ::*UtI3!  *?  J 


Figure  1 1d.  Partial  output  of  BOTSCB  for  problem  1c. 


2.6.4  Problem  2 


Consider  an  asymmetric  wing  section  of  2.76  X  length,  depicted  subse¬ 
quently  in  Figure  17,  with  coordinates  specified  in  Figure  12.  Use  NP  -  17 
points  to  describe  the  body  and  let  NMODE  ■  4.  Compute  the  (^-polarized 
radiation  patterns  for  a  t -polarized  axial  slot  antenna  of  2.06  X  length  and 
0.125  X  width,  centered  about  Z  ■  0,  X  ■  28.796,  and  Y  ■  26.085.  (The  loca¬ 
tion  of  this  slot  corresponds  to  the  eighth  triangle  function  and  is  denoted 
as  antenna  A  in  Figure  17.)  Carry  out  the  analysis  at  10  MHz. 

Solution  -  The  coordinates  for  the  wing  section  are  input  by  the  user.  Then 
the  following-  sequence  of  codes  are  run:  BOTZSS,  BOTINV,  and  BOTRA.  The 
input  data  for  execution  of  these  programs  are  given  in  Figure  12.  The  par¬ 
tial  outputs  from  these  programs  are  shown  in  Figures  13-15.  The  power 
radiation  patterns  are  summarized  in  Figure  16.  These  results  are  plotted  in 
Figure  17  for  antenna  location  A,  as  well  as  locations  B  and  C.  (The  latter 
two  patterns  are  obtained  if  the  aperture  coincides  with  the  first  and  sixth 
triangle  functions  in  Figure  13,  respectively.) 
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Figure  12.  Input  data  for  execution  of  BOTZSS,  BOTINV,  and  BOTRA  for  problem  2. 
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Figure  16a.  Partial  output  of  BOTRA  for  problam  2. 
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Figure  16a.  Partial  output  of  BOTRA  for  problam  2.  (condudad) 


Figure  16b.  Partial  output  of  BOTRA  for  problem  2. 
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Figure  16b.  Partial  output  of  BOTRA  for  problem  2.  ( concluded) 
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Figure  17.  Computed  radiation  patterns  for  wing  section. 


2.6.5  Problem  3 

Consider  a  square  cylinder  with  sides  0.3183  X  and  a  length  of  2.76  X. 

Illuminate  the  cylinder  broadside  (0  “  90°)  with  a  TM  (0-polarized)  wave. 

Repeat  the  analysis  for  a  TE  (ij>-polarized)  wave. 

a)  Compute  the  resulting  bistatic  cross  sections  for  all  polarizations 
both  in  the  horizontal  ($  ■  0°)  plane  and  the  roll  (0  ■  90  )  plane 
at  46 (NT)  equi-spaced  angles. 

b)  Compute  the  magnitude  and  phase  of  the  currents  induced  by  the  TM 
and  TE  illumination. 


i  * 
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Solution  -  The  calculation*  vara  carried  out  at  10  MHz  (X  ■  30  m).  The  body 
was  represented  by  17 (IF)  coordinate  points,  and  4(MM0DE)  modes  were  used. 

The  coordinate  points  are  such  that  the  corners  of  the  cylinder  correspond  to 
triangle  function  peaks.  This  rastriction  is  necessary  to  obtain  proper 
accuracy  in  the  calculations.  The  execution  of  BOTSEG  is  shown  in  Figure  18a 
in  the  time-share  node.  Input  data  to  execute  BOTZSS,  BOTINV,  and  BOTSCB  are 
summarized  in  Figure  18b.  Partial  listings  of  the  inputs  and  outputs  for 
these  programs  are  given  in  Figures  19-21. 

The  results  for  the  two  parts  of  this  problem  are  given  in  Figures  21b 
and  21c.  Comparison  of  the  bistatic  scattering  results  using  the  present 
analysis  with  published  data  is  summarized  in  Figure  IS  of  Volume  I.  The 
aperture  penetration  and  near-field  analysis  options  were  not  exercised  in 
this  problem. 
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Figure  19.  Partial  output  of  BOTZSS  for  problem  3. 
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Figure  21b.  Partial  output  of  BOTSCB  for  problem  3.  (continued) 
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Figure  21b.  Partial  output  of  BOTSCB  for  problem  3.  (concluded) 
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Figure  21c.  Partial  output  of  BOTSCB  for  problem  3.  (concluded) 


2.6.6  Problem  4 

Determine  the  field  penetration  Into  e  fuselage  equipment  compartment 
represented  by  a  BOT,  depicted  In  Figure  22.  Assume  that  the  aperture,  cen 
tered  about  Z  ■  0,  la  1.6S6  X  In  axial  length  and  0.63  X  in  width.  Let  the 
Illumination  of  the  body  be  broadside  (6  -  90°)  at  75  MHz  in  either  8  or  $ 
polarization.  Compute  the  fields  at  the  following  points: 
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4*5.0 

Figure  22.  Equipment  compartment  with  aperture. 


Solution  -  The  solution  was  implemented  by  running  BOTZSS,  BOTINV ,  and  BOTSCB 
in  sequence.  The  parameters  NP  and  NMODE  were  33  and  A,  respectively.  The 
necessary  input  data  to  execute  the  above  programs  are  given  in  Figure  23, 
with  the  partial  outputs  from  BOTZSS,  BOTINV,  and  BOTSCB  shown  in  Figures  24, 
25,  and  26,  respectively.  The  bistatic  scattering  cross  sections  and  the 
currents  on  the  body  with  the  aperture  sealed  are  shown  for  6  and  $  illumina¬ 
tion  in  Figures  26b  and  26c,  respectively.  The  electric  and  magnetic  field 
components  Inside  the  body  at  the  various  sampling  points  are  summarised  in 
Figure  26d. 
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Figure  23.  Input  data  for  execution  of  BOTZSS,  BOTINV,  and  BOTSCB  for  problem  4. 
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Figure  24.  Partial  output  of  BOTZSS  for  problam  4. 
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Figure  24.  Partial  output  of  B0TZSS  for  problem  4.  ( concluded) 
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Figure  26a.  Partial  output  of  BOTSC8  for  problem  4. 
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Figure  26b.  Partial  output  of  BOTSCB  for  problem  4. 


Figure  26b.  Partial  output  of  B0TSC8  for  probtan  4.  (continued) 
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Figure  26b.  Partial  output  of  BOTSCB  for  problem  4.  (concluded) 
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Figura  28c.  Partial  output  of  B0TSC8  for  probtom  4. 
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Figure  26c.  Partial  output  of  BOTSCB  for  problem  4.  (conducted) 
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.  SYSTEMS  SECTION:  DETAILED  PROGRAM  DESCRIPTIONS 


A  detailed  description  of  each  main-line  MM/BOT  program  (Figure  1)  will 
be  given,  including  a  description  of  the  flow  diagram,  subroutine  input/out¬ 
put  arguments,  and  the  special  matrix  properties  used.  A  description  of  the 
common  variables  used  in  the  program,  along  with  storage  methods  used  in 
certain  arrays,  is  given  in  Appendix  A.  A  summary  of  subroutine  calling 
programs  is  contained  in  Appendix  B. 

3.1  BOTZSS  Program 

BOTZSS  generates  the  impedance  submatrices  Z  for  modes  m,n  where  m  ■  0 

mn 

to  NMODE  -  1,  n  •  -m  to  -hn,  and  |m  -  n|  <  NBAND.  The  impedance  matrices  are 

generated  in  the  lower  triangular  portion  of  Z  _.  Symmetry  conditions  are 

then  applied  in  BOTINV  in  order  to  fill  the  entire  Z„_„  matrix.  The  struc- 

oU  1 

ture  of  the  ZBQT  matrix  is  as  follows: 


-1,-1 


Z-1.0 


Z-1,1 


QP7»04I1-107 


where  each  of  the  Z  matrices  is  comprised  of  four  submatrices  as  follows: 
m,n 


V  AD-A087  403 


UNCLASSIFIED 


MCDONNELL  DOUGLAS  RESEARCH  LABS  ST  LOUIS  MO  F/G  20/14 

radiation  and  scattering  from  bodies  of  translation,  volume  II.— ETC (U) 
APR  80  L  N  MEDGYESI-MITSCHANG*  J  M  PUTNAM  F30602-77-C-0233 

RADC-TR-GO-142-VOL-2  NL 


Figure  27  shows  the  flow  diagram  for  BOTZSS.  The  equation  numbers  -refer 

to  the  theoretical  expressions  in  Volume  I, 

The  computation  of  the  Green's  function  kernel  takes  advantage  of  the 

fact  that  G_  is  symmetric  (i.e.,  G  ■  G_  ),  where  t  indicates  the  trans- 
m,n  m,n  t  m,n 

pose  operation.  Only  the  upper  triangular  portion  is  stored  as  indicated  in 

the  appendix;  hence,  (G  ) .  .  is  stored  in  location 

m,n  i,J 

G(i  +  (j-1) j/2)  when  i  <  j 


G(j  +  (i-l)i/2)  when  i  >  j. 


The  section  that  computes  the  impedance  matrix  Z  uses  the  following 

m9n 


symmetries: 


Ztfc  -  (Zfct) 
m9n  t  m9n 


t<Zm!n> 


-  n  Zzt  , 
m,n 


Zzz  -  t(ZZZ  )  . 
m,n  t'  m,n 


Thus,  only  the  upper  triangular  portion  of  each  of  the  ZtC  ,  Zzt  ,  Ztz  ,  and 

n9il  1U)Q  BjQ 

Zzz  needs  to  be  computed.  The  remaining  portion  is  filled  using  the 
m9n 

symmetry  conditions  above. 


In  Figure  27,  (Xp,  Y^)  define  a  point  on  the  BOT,  with  being  the  angle 
subtended  by  the  x-axls  and  the  BOT  generating  curve.  The  parameter  6  -  T/L, 
where  L  is  the  half-length  of  the  BOT  measured  along  the  z-axls  and  T  is  the 
distance  between  adjacent  (Xp,  Yp)  points  on  the  BOT  generating  curve.  Note  a 
triangle  function  subtends  2  T. 


Figure  27.  BOTZSS  flow  diagram. 


3.2  BOTZSS  Subroutines 


3.2.1  Subroutine  CSIMP 

Subroutine  CSIMP  is  e  Slop son  Integration  routine  with  a  calling  state¬ 
ment 

CALL  CS1MP(F,A,B,DEL,IMAX,SI1,S,N,IER). 


/B 

F(x)dx  using  the  method  of  successive  bisections 
A 

of  the  interval  until  either  a  relative  error  of  DEL  is  achieved  or  IMAX 
bisections  have  been  performed.  F  must  be  declared  external  in  the  calling 
program.  The  following  are  returned  by  CSIMP: 

S  -  Approximate  value  of  the  integral. 

SI1  -  Previous  approximation  to  the  integral.  Convergence  has  occurred 


if 


<  DEL. 


N  Number  of  intervals  used  in  computing  S. 

IER  -  Error  return.  IER  «  0  indicates  that  convergence  has  occurred. 


3.2.2  Subroutine  PLOTB 

Subroutine  PLOTB  plots  the  points  on  the  generating  curve  of  the  BOT. 
Points  on  the  BOT  corresponding  to  triangle  function  peaks  are  Indicated  with 
a  plus  sign.  The  calling  statement  is 

CALL  PLOTB(X,Y,N,NR), 

where 

X  -  Array  of  x  coordinates  to  be  plotted 
Y  -  Array  of  y  coordinates  to  be  plotted 
N  -  Number  of  points  to  be  plotted 

NR  -  Number  of  line  printer  rows  to  use  for  the  y-axls. 

The  routine  uses  51  columns  for  the  x-axls,  with  the  dynamic  range  on  both 
the  x  and  y  axes  equal.  Hence,  depending  upon  the  type  of  line  printer  used, 
NR  may  have  to  be  adjusted  in  order  to  obtain  a  plot  that  is  not  distorted 
(i.e.,  the  x  and  y  axes  have  approximately  the  same  physical  length  on  the 
line  printer  output). 


3.3  BOTINV  Program 

BOTINV  fills  the  ZBQT  matrix  using  the  output  file  from  BOTZSS  and 
Inverts  this  matrix  according  to  the  user's  specifications.  The  following 


Figure  28.  ZBqT  matrix  tymmstrtai. 
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Figure  29  shows  the  flow  diagram  for  BOTINV.  Three  types  of  matrix 
Inversions  are  allowed  in  BOTINV:  total  inversion,  main  diagonal  inversion, 
and  partial  inversion.  Each  of  these  options  is  described  next. 


Total  Inversion  -  Total  inversion  is  performed  when  NBAND  >_  2*NM0DE-1.  In 


this  case,  the  matrix  is  stored  by  columns  as  follows: 


(z“) 

m  pti 


i.j 


Z((n  +  NMODE-1) *LS* (2*NM0DE-1) 
+  (j-l)*(2*NMDDE-l)*LS 
+  (m  +  NMODE-1) *LS  +  i) 


(ZZt  )  is  stored  at  ( ZZt  )  +  NM 

i*J  “*  i.J 


(Ztz  )  is  stored  at  (Ztfc  )  +  NM*LS*(2*NM0DE-1) 

m,n  »,n  itj 


.tz 


(Z  )  is  stored  at  (Z„  _)  +  NM 

m.n  m,n  ±  i 
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Once  the  Zg^  matrix  is  filled,  it  is  inverted  using  Gaussian  elimination 
with  partial  pivoting,  and  written  to  disk  file  by  subtutrices. 

Main  Diagonal  Inversion  -  Main  diagonal  inversion  is  performed  when  NBAND  -  L 
In  this  case,  the  individual  diagonal  submatrices  are  inverted  separately 
using  Gaussian  elimination.  For  m  i  0,  the  following  symmetry  is  used: 


ztc 

ZtZ 

ztc 

-  z“ 

— m,  — in 

«_  s. 

-  E»E  _ 

_ a.® 

z*t 

ZZZ 

-z“ 

Zzz 

Hn*-n 

m9m 

which  implies  that 


has  the  same  symmetries.  Thus,  at  a  given  time,  only  Z  where  m  -  0  to 

UyU 

NMODE-1  is  Inverted.  The  symmetries  are  used,  and  the  resulting  (2*NM0DE-1). 
submatrices  are  written  to  disk  file. 

Partial  Inversion  -  Partial  inversion  is  performed  when  NBAND<2*NMODE-l ,  and 
NBAND^l.  In  this  case,  the  ZB(yr  matrix  is  filled  only  with  the  Z^  sub- 
ma trice  for  which  |m-n|  <  NBAND.  The  resulting  ZB()T  matrix  has  a  staircase- 
type  structure.  The  rest  of  the  matrix  is  sparse.  If  each  of  the  Zm  Q 
submatrices  are  thought  of  as  individual  elements,  ZBQT  can  be  considered  as 
a  banded  matrix.  A  modified  LU  decomposition  can  then  be  used  with  all 
arithmetic  operations  replaced  by  the  corresponding  matrix  operations.  The 
result  is  an  L  and  U  matrix  which  are  also  of  a  staircase-type,  but  are  lower 
and  upper  triangular,  respectively,  when  the  submatrices  are  considered  as 
individual  elements.  The  inverse  of  can  then  be  computed  using  forward 

and  backward  substitution,  again  replacing  arithmetic  operations  with  matrix 
operations.  The  result  is  a  full  inverted  ZBq^  matrix  which  is  written  to 
disk  file  by  submatrices. 


The  ZB0T  matrix  is  stored  by  columns ,  if  the  individual  submatrices  Z 
are  considered  as  elements.  Only  the  banded  portion  is  stored.  When  HMODE  “ 
4  and  NBAND  ■  2  (refer  to  Figure  28) ,  ZBQT  is  stored  in  the  following  order: 

Z-3,-3*  Z-2,-3*  Z-3,-2*  Z-2,-2*  Z-l,-2*  Z-2,-l  *  *  * 

2 

Each  submatrix  is  stored  in  LS  successive  locations  by  columns.  For  the 

2 

example  above,  z  2  -3  UDU^^  start  at  index  LS  +1. 

3.4  BOTINV  Subroutines 
3.4.1  Subroutine  LINEQ 

Subroutine  LINEQ  is  a  standard  matrix  inversion  routine  using  Gaussian 
elimination  with  partial  pivoting,  with  the  following  calling  statement  and 
arguments : 

CALL  LINEQ(LL,C,LR) , 

where 

LL  -  Order  of  matrix  to  be  inverted. 

C  -  Array  containing  the  matrix  to  be  inverted,  stored  by  columns.  On 
output,  C  contains  the  inverted  matrix. 

LR  -  Array  of  length  LL  used  as  a  work  space  during  the  pivoting  pro¬ 
cess. 


3.4.2  Subroutine  LIST 

Subroutine  LIST  prints  individual  Y  _  submatrices  on  the  line  printer 

m,n 

and  writes  them  to  a  disk  file.  The  calling  statement  and  arguments  follow: 

CALL  LIST(M,N,Z), 

where 

M  -  Index  m- 
N  -  Index  n 

Z  -  Array  containing  the  Y  submatrix,  stored  by  columns. 

ID|U 


3.4.3  Subroutine  INVBAN 


Subroutine  INVBAN  is  a  modification  of  a  standard  banded  matrix  inver¬ 
sion  routine  using  LU  decomposition  without  pivoting,  where  only  the  banded 

portion  is  stored  by  columns.  Arithmetic  operations  were  replaced  by  their 

2 

corresponding  matrix  operations,  and  indices  were  multiplied  by  LS  since  the 
elements  were  replaced  by  matrices.  The  calling  statement  and  arguments 
follow: 

CALL  INVBAN (LS,NMODE,NBAND,NZ,A,Z, WORK)  , 

where  LS,  NMODE,  and  NBAND  are  described  in  Appendix  A. 

INPUT:  NZ  -  Array  used  for  indexing.  In  a  normal  banded  matrix  A, 

NZ(I)  »  NZ(I-l)  +  (the  number  of  zeroes  below  the  band  in 

column  1-1)  +  (the  number  of  zeroes  above  the  band  in 

column  I),  where  NZ(1)  >0.  If  A  has  order  n,  with  only 

the  banded  portion  stored  by  columns,  then  will  be 

stored  in  location  n(j-l)  +  i  -  NZ(j). 

A  -  Array  containing  the  staircase-type  matrix  to  be  Inverted 

with  storage  details  described  above. 

2 

Z  -  Array  of  length  LS  used  as  work  area. 

WORK  -  Array  of  length  LS  used  as  work  area. 

In  addition,  three  variables  are  passed  in  common  as  follows: 

COMMON  NM,JK(4) ,LR 

where 

NM  -  Number  of  triangle  functions 
JK  -  Work  array  of  length  4 
LR  -  Work  array  of  length  LS. 

3.4.4  Additional  Subroutines 

The  subroutines  MULTS,  MULT,  REPLACE,  and  ZERO  perform  matrix  operations 
which  are  documented  in  the  listings. 


3.5  BOTRA  Program 

BOTRA  computes  the  radiated  far  and  near  fields  produced  by  one  or  more 
slot  antennas  on  the  BOT  (see  Section  5  of  Volume  I).  Figure  30  shows  the 
flow  diagram  for  BOTRA  with  equation  numbers  referenced  to  Volume  I. 

The  parameters  (p^,  <j>p)  define  (X^,  Y^)  in  cylindrical  coordinates. 
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Figure  30.  BOTRA  flow  diagram. 
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Figure  30.  BOTRA  Row  diagram.  ( conti nuad) 
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Figure  30.  BOTRA  flow  diagram,  (concluded) 


3 . 6  BOTRA  Subroutines 


3.6.1  Subroutine  PLANE 

Subroutine  PLANE  computes  the  mode  Independent  part  of  the  transfer 
matrices  given  in  Equations  (34-37)  of  Volume  I.  The  calling  statement  and 
argument  list  is  given  below: 

CALL  PLANE ( RT ,RP ,TH, NT, PHI) , 

where 

RT  -  Array  containing  the  mode  independent  part  of  the  8  polarized 
transfer  matrices,  for  NT  different  angles.  See  variable  RT  in 
Appendix  A  for  storage  details. 

RP  -  Array  containing  the  mode  independent  part  of  the  <f>  polarized 
transfer  matrices,  for  NT  different  angles.  See  variable  RP  in 
Appendix  A  for  storage  details. 

INPUT:  TH  -  Array  containing  the  NT  8  angles. 

NT  -  Number  of  transfer  matrices  to  be  calculated,  corresponding 
to  the  angles  in  arrays  TH  and  PHI. 

PHI  -  Array  containing  the  NT  <p  angles. 

In  addition,  the  body  parameters  are  passed  to  subroutine  PLANE  in  the  common 
block  BODY.  The  parameters  in  BODY  are  described  in  Appendix  A. 

3.6.2  Subroutine  NEARB 

Subroutine  NEARB  calculates  the  electric  and  magnetic  modal  current 
coefficients  for  the  near-field  calculations  (see  Section  6  of  Volume  I). 

The  flow  diagram  for  NEARB  is  shown  in  Figure  31.  The  calling  statement  and 
argument  list  is  described  below: 

CALL  NEARB (XTEST, YTEST, ZTEST, ZM) , 

where 

INPUT:  XTEST  -  x  coordinate  of  test  point. 

YTEST  -  y  coordinate  of  test  point. 

ZTEST  -  z  coordinate  of  test  point. 

OUTPUT:  ZM  -  Array  of  modal  current  coefficients  for  mode  M.  Storage 
details  are  given  in  the  appendix. 


In  addition,  other  variables  are  passed  to  NEARB  in  common  blocks  BODY,  TEST, 
and  INT.  These  variables  are  described  in  Appendix  A. 


Figure  31 .  Subroutine  NE ARB  flow  diagram. 


3.6.3  Subroutine  LPCUR 

Subroutine  LPCUR  lists  and  plots  the  BOT  currents.  The  entire  array  of 
t  and  z  directed  nodal  currents  is  first  printed,  followed  by  plots  of  the  t 
and  z  directed  currents  (magnitude  and  phase  as  a  function  of  the  normalized 
z  coordinate)  for  each  triangle  function  on  the  body.  The  calling  statement 
and  argument  list  follows: 

CALL  LPCUR (NM0DE,NM, CUR) , 

where  all  of  the  input  variables  are  described  in  Appendix  A. 

3.6.4  Subroutine  PLOT 

Subroutine  PLOT  plots  the  magnitude  and  phase  of  the  currents  on  a  given 
triangle  function.  The  calling  statement  is  as  follows: 

CALL  PL0T(Y1,Y2),  . 


"SWWIWWSwse  fc-S***  O  .  ■< 


where 

INPUT:  Y1  -  Array  containing  the  current  magnitude  at  41  equally  spaced 
z  coordinates. 

Y2  -  Array  containing  the  current  phaae  at  41  equally  spaced  z 
coordinates. 


3.7  BOTSCM  Program 

BOTSCM  computes  scattered  far  fields  in  the  monostatic  mode  (see  Section 
5.2  of  Volpme  I).  Figure  32  shows  the  flow  diagram  for  BOTSCM  with  equation 
numbers  referenced  to  Volume  I.  All  of  the  subroutines  called  by  BOTSCM  are 
described  in  the  BOTRA  subroutine  section.  In  the  flow  diagram,  all  transfer 
matrices  are  computed  in  one  call  to  PLANE,  while  in  BOTRA  and  BOTSCM,  one 
transfer  matrix  is  calculated  at  a  time. 


3.8  BOTSCB  Program 

BOTSCB  computes  scattered  far  fields  in  the  bistatlc  mode,  after  which 
scattered  near  fields  can  be  calculated  with  or  without  an  aperture  present. 
Figure  33  shows  the  flow  diagram  for  BOTSCB  with  equation  numbers  referenced 
to  Volume  I. 


■ 
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3.9  BOTSCB  Subroutines 

Subroutines  PLANE,  NEARB,  LPCUR,  and  PLOT  are  described  in  the  BOTRA 
subroutine  section.  Subroutine  LINEQ  is  described  in  the  BOTINV  subroutine 
section. 

3.9.1  Subroutine  APPAR 

Subroutine  APPAR  calculates  the  aperture  parameters  LT,  VLOW,  VHGH,  IT, 
and  NAM  which  are  described  in  Appendix  A.  These  variables  are  returned  in 
common  block  S.  The  calling  statement  and  input  parameter  list  is  as  follows: 

CALL  APPAR(X0,Y0,X1,Y1,XH,YH,DH,NP, IEDGE) , 

where 

XO,  YO,  XI,  Y1  -  Described  in  the  user  input  section  (Section  2). 

XH,  YH,  DH,  NP,  IEDGE  -  Described  in  Appendix  A. 
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Figure  32.  BOTSCM  flow  diagram . 
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(Equation  (61)) 


(Equation  (51)] 
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Figure  32.  BOTSCM  flow  diagram,  (condudad) 


CD 

Wfura  S3.  BOT8C8  flow  diagram.  ( continued) 


Figure  33.  B0TSC8  flow  diagram,  (concluded) 

3.9.2  Subroutine  APULSE 

Subroutine  APULSE  computes  the  t-  and  t-dlrected  aperture  pul.ee  func¬ 
tions.  The  calling  statement  and  argument  list  Is  described  below: 

CALL  APPAR(FPP(1, J) ,FPPZ(1, J) ,IT(J) ,LT(J) ,VLOW(J) ,VHGH(J) ,DH) , 

where 

OUTPUT:  FPP(1, J) ,PPP(2, J) ,FPP(3, J)  -  The  three  t-dlrected  aperture 

pulse  functions  corresponding  to  aperture  triangle  func¬ 
tion  J  are  returned  in  these  three  successive  storage 


locations.  Saa  variable  FPP  in  tha  appendix  for  storage 
details. 

FPPZ(1,J) ,FPPZ(2, J) ,FPPZ(3,J)  -  The  three  z-dlrected  aperture 
pulae  functions  corresponding  to  aperture  triangle  func¬ 
tion  J  are  returned  in  these  three  successive  storage 
locations. 

INPUT:  IT(J) ,LT(J) ,VLOW(J) ,VHGH(J)  -  Aperture  parameters  corresponding 

to  the  J-th  aperture  triangle  function,  as  described  in 
Appendix  A. 

DH  -  Array  containing  BOT  segment  lengths. 

3.9.3  Subroutine  ASYMAC 

Subroutine  ASYMAC  computes  the  equivalent  BOT  currents  in  the  presence 
of  an  asymmetric  aperture  (the  underlying  analysis  is  given  in  Section  7  of 
Volume  I).  The  flow  diagram  for  ASYMAC  is  shown  in  Figure  34.  The  calling 
statement  and  argument  list  is  described  below: 

CALL  ASYMAC (NMDDE,Z0,Z1,Y, CURT, CURP, VT, VP ,NM,DH) , 

where 

INPUT:  NMODE,  ZO,  Z1  -  Described  in  user  input  section. 
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Y  -  Work  array  of  length  LS  used  to  store  the  Y  sub- 

m,n 

matrices  as  they  are  read  from  disk  file. 

CURT  -  Array  of  modal  currents  (resulting  from  a  0-polarised 
Incident  wave)  on  the  BOT  without  an  aperture  present. 

On  returning  from  ASYMAC,  the  array  CURT  will  contain  the 
modal  currents  on  the  BOT  with  the  aperture  present. 

CURP  -  Array  of  modal  currents  (resulting  from  a  ^-polarized 
incident  wave)  on  the  BOT  without  an  aperture  present. 

On  returning  from  ASYMAC,  the  array  CURP  will  contain  the 
modal  currents  on  the  BOT  with  the  aperture  present. 

VT  -  Work  array  of  length  LS. 

VP  -  Work  array  of  length  LS. 

NM,  DH  -  Described  in  Appendix  A. 


APPENDIX  A:  DICTIONARY  OF  COMMON  PROGRAM  VARIABLES 
(Input  variables  are  described  in  Section  2.3;  equation  numbers  refer 
to  expressions  in  Volume  1;  page  numbers  refer  to  Volume  II) 


ANG 

BK 

BKL 

CAP 


CAT 


CP  (I) 

CUR 


CURP 


CURT 


CV(I) 


DELTA 

DH(I) 

DIOR 


Input  array  of  fixed  radiation  or  scattering  angles  (p.  8). 

Wave  number  (meters-*) . 

BK*L. 

Array  containing  equivalent  aperture  currents  resulting  from  a 

^-polarized  incident  wave  [Equation  (86)].  CAP(J)  contains  the 

t-dlrected  current  on  the  J-th  aperture  triangle  function  (i.e. , 

C  ).  The  z-directed  current  (i.e.,  C  )  is  stored  in  CAP(J  + 
ap  ap 

NAM). 

Array  containing  equivalent  aperture  currents  resulting  from  a 
0-polarized  incident  wave  [Equation  (86) ] .  See  variable  CAP  for 
storage  details. 

Cosine  of  4>  angle  for  generating  curve  segment  I.  (See  Figure  5 
of  Volume  I.)  Corresponds  to  cos$q  in  Equations  (38)-(41). 

Array  containing  modal  t-  and  z-dlrected  currents  on  the  BOT 
[Equation  (10)].  CUR  [(m  +  NMODE  -  1)*LS  +  J]  contains  the  t- 
directed  current  for  mode  m  on  triangle  function  J.  The  z- 
directed  current  is  stored  in  CUR  [ (m  +  NMODE  -  1) *LS  +  J  +  NM] . 
Used  in  BOTRA. 

Array  containing  modal  t-  and  z-directed  currents  on  the  BOT  for  a 
^-polarized  incident  wave  [Equation  (10) ] .  See  variable  CUR  for 
storage  details.  Used  in  BOTSCB. 

Array  containing  modal  t-  and  z-dlrected  currents  on  the  BOT  for  a 
0-polarized  incident  wave  [Equation  (10)].  See  variable  CUR  for 
storage  details.  Used  in  BOTSCB. 

Cosine  of  angle  for  generating  curve  segment  I.  (See  Figure  5 
of  Volume  I.)  [Equation  (17)].  CV(J)  corresponds  to  on 
segment  J. 

GAMMA/L. 

Length  of  generating  curve  segment  I  (meters) . 

x/180#. 


ESC 


ESCP 

ESCT 

EO 

ETA 

EVP 

EVT 

FPP 

FPPZ 

G 

GAMMA 

GP 

GT 


-  Array  containing  electric  near -fie Id  radiation  components  of  £(r') 
in  the  x,  y,  and  z  directions  [Equation  (67)],  stored  in  ESC (1-3) , 
respectively.  Used  in  BOTRA. 

-  Array  containing  electric  near-field  scattering  components  result¬ 
ing  from  a  ^-polarized  incident  wave  [Equation  (67)].  Used  in 
BOTSCB . 

-  Array  containing  electric  near-field  scattering  components  result¬ 
ing  from  a  6-polarized  incident  wave  [Equation  (67)].  Used  in 
BOTSCB. 

-  Input  array  for  slot  antenna  excitation  [Equation  (46)].  (p.  9) 

-  n  -  /y/e  -  376.707  fl. 

ro  o 

-  Array  containing  aperture  E-flelds  resulting  from  a  ^-polarized 

incident  wave  [Equation  (89)].  EVP(J)  contains  the  t-dlrected 

E-field  on  the  J-th  aperture  triangle  function.  The  z-dlrected 

E-field  is  stored  in  EVP(J  +  NAM).  EVP(J)  -  EV*;  EVP(J  +  NAM)  - 
~  q 
EVZ. 

q 

-  Array  containing  aperture  E-flelds  resulting  from  a  8-polarized 
incident  wave  [Equation  (89) ] .  See  variable  EVP  for  storage 
details. 

-  Matrix  containing  the  t-directed  pulse  functions  on  the  aperture, 

F®  with  a  -  t  [Equation  (80)].  FPP (I  »  1  to  3,  J)  contains  the 
3  *  q  .  t  ^ 

three  pulse  functions  F.  ,  F,  ,  and  F.  .  ,  for  the  J-th 

3  _  IfP  J *P  J  1*P 

aperture  triangle  function. 

-  Matrix  containing  the  z-dlrected  pulse  functions  on  the  aperture, 
F?  with  a  *  z  [Equation  (80) ] .  See  variable  FPP  for  storage 

j  »q 

details . 

-  Array  contalnihg  the  integrated  Green's  function  kernel  G 

u»n 

[Equation  (21)].  In  BOTZSS,  G  is  symmetric  with  only  the  upper 
triangular  portion  stored  by  columns  from  index  1  to  (NP  -  1)* 
NP/2.  G.  ,  is  stored  in  location  G[i  +  (J  -  l)*J/2]  when  i  <_  J. 

J 

-  Half  length  of  a  triangle  function  base  (meters) . 

-  Array  containing  the  fields  for  the  NT  radiation  angles  [Equetlon 
(43)]  to  compute  ^-polarized  gain. 

-  Array  containing  the  fields  for  the  NT  radiation  angles  [Equetlon 
(43)]  to  compute  9-polarlzed  gain. 


HO 


-  Array  containing  the  integrated  Green's  function  kernel  HQ(q) , 

(Y  **  0)  [Equation  (63)]  for  the  magnetic  near  fields  [Equations 
(63)-(66)]. 

V 

HI  -  Array  containing  the  Integrated  Green's  function  kernel  H^Cq) , 

(Y  **  1)  [Equation  (63)]  for  the  magnetic  near  fields  [Equations 
(63) —(66) ] . 

HSC  -  Array  containing  magnetic  near -field  radiation  components  in  the 

x,  y,  and  z  directions  [Equation  (71)].  HSC  ■  H(r'),  used  in 
BOTRA. 

HSCP  -  Array  containing  magnetic  near-field  scattering  components  result 
ing  from  a  ^-polarized  incident  wave  [Equation  (71)].  HSCP  ■ 
H(r'),  used  in  BOTSCB. 

HSCT  -  Array  containing  magnetic  near-field  scattering  components  result 
ing  from  a  9-polarlzed  incident  wave  [Equation  (71)].  HSCT  ■ 
H(r'),  used  in  BOTSCB. 

1EDGE  -  Indicates  whether  the  generating  curve  is  open  or  closed.  IEDGE 
0  for  closed  and  1  for  open. 

IPLANE  -  Input  array  indicating  whether  corresponding  element  of  array  ANG 
is  a  0  or  <j>  angle,  (p.  8). 

IS  -  Input  array  for  specifying  location  of  slot  antennas.  [See 

Equations  (46)— (49) . ]  (p.  9). 

IT(K)  -  Indicates  whether  aperture  segment  K  starts  or  terminates  the 
aperture.  IT  -  -  1  for  the  start,  IT  -  0  for  interior,  and 
IT  ■  +  1  for  termination. 

KG  -  NP  -  1. 

L  -  Half  length  of  the  BOT  (meters) . 

LS  -  Order  of  each  Z  submatrix.  LS  »  NP  -  3. 

m,n 

LSS  -  LS*LS 

LT(K)  -  Triangle  function  peak  number  for  aperture  segment  K. 

M  -  Mode  number  m. 

MC  -  Input  variable.  M  in  Equation  (26).  (p.  7). 

N  -  Mode  number  n. 

NAM  -  Number  of  triangle  function  peaks  on  the  aperture  portion  of  the 
generating  curve.  NAM  should  be  greater  than  1,  if  an  aperture 
is  present. 


NAM2  -  NAM*2. 

NANG  -  Input  variable.  Number  of  fixed  radiation  or  scattering  angles, 

(p.  8). 

NBAND  -  Input  variable.  Number  of  submatrix  diagonal  bands  used  in 

(P-  7). 

NM  -  Number  of  triangle  functions  on  the  BOT  generating  curve.  NM  “ 

(NP  -  3)/2. 

NMODE  -  Input  variable.  Number  of  non-negative  modes,  (p.  6). 

NM2  -  Order  of  each  Z  submatrix.  NM2  ■  NP  -  3. 

m,n 

NM4  -  NM*4. 

NP  -  Number  of  points,  on  the  BOT  generating  curve.  (See  Section  2.4.) 

NPT  -  Input  variable.  Number  of  diagonal  bands  used  in  Impedance 

matrices,  (p.  7). 

NSA  -  Number  of  slot  antennas  on  the  BOT. 

NT  -  Input  variable.  Number  of  radiation  and  scattering  angles. 

(p.  8). 

NTEST  -  Input  variable.  Number  of  test  points  for  near-fields,  (p.  10). 

PHII  -  <j>  angle  for  the  incident  wave  (degrees).  PHII  ■  <p±  in  Equation 

(51) . 

PHIR(K)  -  <p  angle  for  the  radiated  fields  (degrees) .  PHIR(K)  -  <t>r  in 
Equation  (38) . 

PHIS(K)  -  <p  angle  for  the  scattered  fields  (degrees).  PHIS(K)  -  $  in 

s 

Equation  (51) . 

R(I)  -  Distance  from  origin  to  generating  curve  segment  I  (meters) 

[used  in  Ar  in  Equations  (38) -(41)]. 

^  td>  2  4) 

RP  -  Array  containing  the  mode  Independent  portion  of  the  R  and  R 

n  n 

matrices  (l.e.,  with  the  a  term  removed),  resulting  from  a  4>- 
polarized  incident  wave  [Equations  (40)-(41)].  The  order  of 
storage  is 

|(R^)it  i  -  1  to  NM  followed  by 
|(R**)i.  i  -  1  to  NM 

RP  may  contain  the  transfer  matrices  for  several  (♦,  6) 
angles.  In  this  case,  the  starting  index  is  off-set  by  a  multi¬ 
ple  of  2*NM. 


c6  26 

Array  containing  the  mode-independent  portion  of  the  R  and  R 

n  n 

matrices  (i.e.,  with  the  a  term  removed),  resulting  from  a  0- 
polarized  incident  wave  [Equations  (38) —(39) ] .  The  order  of 


1  to  NM 


storage  is 

{(Rn0)i*  1  *  1  t0  followed  by 
{(Rn9)i*  1  "  1  to  m) 

As  in  RP,  the  starting  index  may  be  off-set  by  a  multiple 
of  2*NM. 

SP(I)  -  Sine  of  <p  angle  for  generating  curve  segment  .1.  (See  Figure  5  of 

Volume  I.)  Corresponds  to  sin<t>  in  Equations  (38),  (40). 

11  <1 

SPP  -  Array  containing  0W  for  the  NT  scattering  angles  [Equation  (51) ] . 

SPT  -  Array  containing  for  the  NT  scattering  angles  [Equation  (51)]. 

STP  -  Array  containing  for  the  NT  scattering  angles  [Equation  (51) ] . 
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STT  -  Array  containing  a  for  the  NT  scattering  angles  [Equation  (51)]. 

SV(I)  -  Sine  of  angle  for  generating  curve  segment  I.  (See  Figure  5  of 

Volume  I.)  [Equation  (17).]  SV(J)  corresponds  to  v  on  segment  J. 

^  £ 

T  -  Array  containing  the  values  of  the  triangle  functions  T^. 

T[(K  -  1)*4  +  p]  contains  the  values  of  the  kth  triangle  function 
over  the  p-th  segment  forming  it.  1  £  P  £  4.  [Equations  (17) — 

(20).]  (p.  6). 

TEXC  -  Input  array  indicates  t-excitation  on  slot  antenna,  (p.  9). 

THI  -  6  angle  for  the  incident  wave  (degrees) .  THI  -  0^  in  Equation 

(51). 

THR(K)  -  0  angle  for  the  radiated  fields  (degrees).  THR(K)  -  0r  in 

Equation  (43) . 

THS(K)  -  6  angle  for  the  scattered  fields  (degrees)  [Equation  (51)]. 

•  £ 

TP  -  Array  containing  the  values  of  T^.  The  storage  method  is  the 

same  as  for  T  [Equations  (17)-(20)].  (p.  6) 

TZ  -  Array  containing  the  values  of  T*.  The  storage  method  is  the 

same  as  for  T  [Equations  (17)-(20)]. 

U  -  Imaginary  number  i. 

UMN  -  Array  containing  values  of  the  u  function  [Equation  (23)]  needed 

mn 

in  numerical  integration  of  the  Green's  function  kernel.  [See 
Equation  (27) . ] 


'%  • 


VHGH(I) 

VLOW(I) 

VM 


VN 

VP 


VT 


XO 

XI 

XH 

XS(I) 

XTEST 

Y 


YO 

Y1 

YA 

YH 

YP 


Upper  "part"  width  asaoclatad  with  I-th  aperture  port. 

Lower  "part"  width  aaaoclated  with  I-th  aperture  port. 

Array  of  voltages  corresponding  to  mode  M.  VM(K)  contains  t- 
directed  voltages  V* .  on  triangle  function  K.  VM(K  +  HM)  contains 
z-directed  voltages  V  ^  on  triangle  function  K  [Equation  (11)). 
Array  of  voltages  corresponding  to  mode  N.  The  storage  method  is 
the  same  as  for  VM  [Equation  (11) ] . 

Array  containing  the  equivalent  BOT  voltages  with  an  aperture 
present,  resulting  from  a  ^-polarized  incident  wave  [Equation 
(90)].  VP(J)  is  VC.  and  contains  the  t-directed  voltage  on  the 

BOT  triangle  function  J.  The  z-directed  voltage  V  is  stored  in 

■3 

VP(J  +  NM) . 

Array  containing  the  equivalent  BOT  voltages  with  an  aperture 
present,  resulting  from  a  6-polarlzed  incident  wave  [Equation 
(90) ] .  VT(J)  contains  the  t-directed  voltage  on  the  BOT  triangle 
function  J.  The  z-directed  voltage  is  stored  in  VT(J  4-  MM). 
Starting  x  coordinate  for  aperture,  (p.  10). 

Ending  x  coordinate  for  aperture,  (p.  10). 

Input  array  of  x  coordinates  for  BOT.  (p.  7). 
x  coordinate  for  generating  curve  segment  I  (meters)  [Equation 
(28)]. 

Input  variable  for  near-field  test  point  x'  in  r'  [Equations  (67), 
(71)].  (p.  10). 

Array  containing  the  Y  submatrix  [Equations  (43) ,  (51) ] .  In 

m,n 

the  near-field  analysis,  however,  Y  contains  the  measurement 

matrix  ZM.  Y  is  stored  by  columns. 
m,n 

Starting  y  coordinate  for  aperture,  (p.  10). 

Ending  y  coordinate  for  aperture,  (p.  10). 

Array  containing  the  aperture  admittance  matrix,  stored  by  columns 
[Equation  (88)].  YA  is  of  order  MAM. 

Input  array  of  y  coordinates  for  BOT.  (p.  7). 

Array  used  as  intermediate  storage,  containing  r£($b,  ♦#)  n 
for  a  given  m,  n,  and  scattering  angle  [Equation  (51)]. 
y  coordinate  for  generating  curve  segment  1  (meters)  [Equation 
(28)]. 


YS(I) 


YT 


YTEST 

Z 

ZEXC 

ZM 


ZTEST 

ZO 

Z1 


g 

-  Array  used  aa  Intermediate  storage,  containing  R($  ,  0  )  Y 

MS  S  tt»a 

for  a  given  m,  n,  and  scattering  angle  [Equation  (51)]. 

-  Input  variable  for  near-field  test  point  y1  in  r'  [Equations  (67), 

(71)].  (p.  10). 

-  Array  containing  the  Z  submatrix  [Equations  (10),  (17)-(20)]. 

m,n 

In  the  program  BOTINV,  however,  Z  contains  the  entire  Z^^, 
matrix.  Z  is  stored  by  columns. 

D|Il 

-  Input  array  [Equation  (48)]. 

-  Artay  containing  the  electric  and  magnetic  modal  current  coeffi¬ 
cients  for  near-field  calculations  [Equations  (67),  (71)].  ZM 
uses  the  same  storage  location  as  the  Y  matrix,  which  is  of  order 
LS,  stored  by  columns.  Rows  1  through  3  of  ZM  contain  the  M-th 
modal  current  coefficients  for  the  electric  near-field  components 
in  the  x,  y,  and  z,  respectively,  at  the  point  (XTEST,  YTEST, 
ZTEST).  Similarly,  rows  4  through  6  of  ZM  contain  the  M-th  modal 
current  coefficients  for  the  magnetic  near-field  components  in 
the  x,  y,  and  z  directions,  respectively. 

-  Input  variable  for  near-field  test  point  z'  in  r'  [Equations 

(67),  (71)].  (p.  10). 

-  Input.  ZO  has  different  definitions  in  the  programs  BOTRA  and 

BOTSCB  [Equation  (49)].  (p.  10). 

-  Input.  Z1  has  different  definitions  in  the  programs  BOTRA  and 

BOTSCB  [Equation  (49)].  (p.  10). 


APPENDIX  B:  SUBROUTINE  CALLING  PROGRAMS  (BOTZSS, 
BOTIN V,  BOTRA,  BOTSCM,  AND  BOTSCB  ARE 
MAIN  PROGRAMS) 


Subroutine 

Calking  program 

APPAR 

BOTSC8 

APULSE 

ASYMAC 

ASYMAC 

BOTSCB 

CSIMP 

BOTZSS 

DBCON 

BOTRA  BOTSCB,  BOTSCM 

FUNC 

CSIMP 

INVBAN 

BOTIN  V 

LINEQ 

ASYMAC,  BOTIN V,  INVBAN 

LIST 

BOTINV.  INVBAN 

LPCUR 

BOTRA  BOTSCB 

MULT 

INVBAN 

MULTS 

INVBAN 

NEARB 

BOTRA  BOTSC8 

PLANE 

BOTRA  BORSC8.  BOTSCM 

PLOT 

LPCUR 

PLOTS 

BOTZSS 

REPLACE 

INVBAN 

SINC 

BOTRA  BOTSCB,  BOTSCM 

ZERO 

INVBAN 

